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Abstract 

We develop a Radon like transformation, in which P quasiprobability 
distribution for spin 1/2 states is written in terms of the tomographic 
probability distribution w. 
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Introduction 

The idea of finding a description of quantum states in terms of (generalized) 
functions which are analogous to probability distributions in particle phase space 
was realized by Wigner Q (see also ref.Q). The Wigner function, which is a real 
valued distribution, turned out to take negative values in some domains of phase 
space for some quantum states. Due to this the Wigner function was called a 
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quasiprobability distribution (or a quasidistribution) . The quantum mechanical 
basic equations were rewritten in terms of the Wigner function by MoyalH. 




FeynmanQ] has even discussed the possibihty to drop the assumption that the 
probabihties of an event must be a nonnegative number. These attempts are 
compared with the permanent wish to explain quantum mechanics in intuitively 
more acceptable classical notions. The hidden variables theories are examples 



Addressing the question of how to make quantum mechanics look like a 
hidden variable theory and vice versa, an analog of Glauber-Sudarshan[^ P 
distribution was constructed for spin 1/2 states by one of the authors 0. The 
introduced quasiprobability distribution was expressed in terms of the conven- 
tional density matrix description of spin 1/2 states. This distribution was ap- 
plied to propose a hidden variable theory, whose predictions agree with quantum 
mechanical ones in various cases (for a tutorial, see ref.[||). This technique was 
generalized to any pair of noncommuting observables and used to reinterpret 
Einstein-Podolsky-Rosen-Bohm "paradox"^, and the role of Feynman's neg- 
ative probabilities was made concrete in the example of interferometers with 
which path (WW for the German welcher Weg) detectors|l^. Furthermore, 
the question of operator ordering was addressed in ref . , where some analo- 
gous distributions were proposed. 

Recently the tomographic methods of measuring quantum states were dis- 
cussed jl^. In these schemes the invertible Radon transform was used to ex- 
press the Wigner function in terms of positive probability distributions (called 
marginal distributions). Such marginal distributions can be obtained directly 
from measurements, and are used to completely characterize a quantum state. 
The tomography of spin states was discussed in refs.||l^. On the other hand, 
for spin 1/2 states the relation between quasiprobability P distribution^ and 
positive probability distribution has not been discussed yet. 

The aim of our work is to find the map which is analog of Radon transform 
used in tomographic schemes jl^ for connecting the quasiprobability distribution 
for spin 1/2 particles with measurable probability distributions. We also express 
the density matrix in terms of the distributions P and w, and relate the last 
with Bloch vector. 

This Letter is organized as follows: in section |l| the properties of P quasiprob- 
ability distribution are reviewed; in section ^ the w probability distribution is 
discussed. In section || their relation is shown. The paper ends with a few 
concluding remarks and some calculational tools are collected in the appendix. 



1 A Review on P (s) Spin 1/2 Quasiprobability 
Distribution 



Some years ago, one of the authors has introduced a quasiprobability distribu- 
tion P (s) in order to describe a spin 1/2 quantum state[Q. In some features 
it can be viewed as an analog of Glauber-Sudarshan P distributions for light 




of attempts in this direction 



2 



(equivalently harmonic oscillator) quantum statesg. For a tutorial, the reader 

is referred to ref. 

citeScuMex. 

In order to define P (s) we should first introduce the S operator given by 



dx 



where O is an operator (e.g. Pauli spin operator ctJ) and c is a real variable 
which plays the role of its classical analog (e.g. s^). Then, given an operator 
Q {a) , we can define a classical counterpart Q (s) by 

Q{a) ^ J (fsQ {sa;,Sy, Sz) S (s^; -a^)S {sy -ay)S {s^ - a^) , 

where should be noticed that an operator ordering was chosen. Quantum mean 
value 

'Q)^Tr(pQ), (1) 



and "classical" mean value 

{Q)=Jd'sP{s)Q{s), (2) 

suggest the definition 

P[s) ^Tr [p5 {sx -ax)S {sy -ay)S (s^ - a^)] , (3) 

which makes quantum and "classical" mean values (eqs. (|^) and (|^)) to coincide. 

For the purposes of ref. (treat Stern-Gerlach like experiments) the referred 
author has "traced out" the y variable. We will not proceed in this way! We 
will treat P (s) in its full glory, in the same spirit as in ref. |pl| , where some 
important properties are made more evident. 

1.1 An Equivalent P (s) Presentation 

With a few algebraic manipulation, we get another expression for P (s), in which 
some important features easily appears. 

Just inserting the definitions of S {rrii — ai) and using eigenvectors of to 
calculate trace, we get 

^ (^^) = E / § / ^ / ^e-'^^'^e-^'^'^e-'^^' (a.| pe^^^-^ e^^^^ e^^^' |a,) , (4) 

a 

where 

a^la^) = a \a^) . 

Standard quantum mechanical calculation immediately leads to 

{a^lpe'^^-e^^^ye'"^' \a:,) = ^ e*'"'e*«''e'>^^ (5) 



b,c 



{C. 



■x I by) (by I a^) (a^ \p\ c^) 
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where 

Ox\Cx) = c\Cx) , 

^\hy) = h\by) . 

Substitution of eq. (||) in eq. (i) yields 

^ (^0 = E («^- b, a)) (c, I by) {by I a,) (a, |p| c,) , 

a.h.c 

where we made use of the notation 

5^ (s - (c, 5, a)) = (5 (sa; - c) (5 {sy - b) 5 {sz - a) . 
Decomposition of \cx) in \az) basis leads to 

^(^1 = E<53(s-(c,6,a))(c. |a.) (6) 

{(a^ I Cx) {az \p\ ttz) + (~az \ c^) {az \p\ - a^)} . 

The above expression (which corresponds to eq.(2.18) in ref.|jri|, but written 
in coordinates) makes clear the singular character of P (s) for quantum mechan- 
ical allowed states. It can assume non null values only in the eight vertices of a 
cube centered at the origin and parallel to the Cartesian axis. This is an obvious 
manifestation of the fact that only eigenvalues are allowed as results of a single 
measurement in quantum mechanics. 

We make use of this singular character and define p (c, b, a) by 

P (sO = E is- (c, b,a))p (c, b, a) . (7) 

a,b,c 

Straightforward calculation (see appendix) leads to 







+ i) [P++ 


+ P+- 






- i) [P++ 


- P+- 


= 




- i) [P++ 


+ P+- 






+ i) [P++ 


- P+- 


= 




-^)[P- 


+ P-+ 






+ i)[p— 


- P-+ 






+ i)[p— 


+ P-+ 






-^)[P- 


- P-+ 
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where we used the simple notation 



P++ = { + z\P\+z), 

p+- = 



P-+ = {-z \P\ 



zl , 



1.2 Density Matrix from P (s) 

We now show how to invert relation and obtain density matrix from P (s) 
quasiprobability distribution. Then, with aid of simple properties of density 
matrices, we derive some properties which a P (s) distribution should obey to 
correspond to a quantum mechanical possible state. 
From the first two lines of eqs. (^) one easily obtains 

P++ = (l-^)p(l,l,l) + (l + ^)p(-l,l,l), (9) 

P+- = + 

and using the properties of hermiticity and unitary trace of density matrix, one 
gets 

P++ P4 
P*+- 1 - P++ 

This result shows us that from an information theoretical viewpoint all other 
six p (c, b, a) are redundant information. Algebraically, this redundant informa- 
tion manifests itself as a set of constraints which a P (s) distribution must obey, 
in order to represent a quantum mechanical allowed state (e.g. one which obeys 
uncertainty relations). One of this constraints is just marginal distribution 
property: 

(c, 6, a) = (a^ \p\ a^) ^Vis^^a), 

b,c 

with similar relations holding for V {s^ = c) and V {sy = 6), where V (*) denotes 
probability of *. So now it is clear that P (s) is a complex valued quasiprobability 
distribution. 



1.3 Examples 

We now work out some simple examples. 
i) Spin up in z-axis; 

As its density matrix in \az) basis is given by 



P+z 



1 
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one immediately obtains 



p{-l,-l,l) = -{l + i), 



p(-l,l,l) = p(l,-l,l) = -(l-i), 
p(c,6,-l) = 0. 

This example was also worked out in ref. Q, but considering only x and z 
directions. If we add all y possibilities (6 = ±1) we reobtain that result. 
ii) Spin "up" in x-axis: 
Equivalently, the density matrix is given by 



P+x 



and then 



P(l,l,l) 

p{-l,b, a) 

in) Spin "up" in y-axis: 
Density matrix is given by 



1/2 1/2 
1/2 1/2 



0. 



-4(1 + ^), 



and then 



P+y = 
p{c,l,a) = 



1/2 

i/2 



1 

4 



-i/2 
1/2 



= p(-l,-l,-l) = — , 

= p(l,-l,-l) = ^. 

Those examples show how peculiar is the behaviour of P (s) quasiprobability 
with respect to cube symmetries. This is a consequence of the noncommutativity 
of S operators in eq. (||), which makes P(s) distribution noncovariant. In case 
we had used the symmetric Margenau-Hill P, Ps{s), or the Wigner-Weyl P, 
Pw (s), as defined in ref . |ll[| , they would have a covariant behaviour, but in 
this work we intend to discuss the original P (s) distribution (denoted Pxyz in 
ref.0). 

iv) Unpolarized spin: 

In this last example, density matrix is given by 



Po 



1/2 
1/2 
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and one gets 

= = i(l-i). 

Those four examples span all possibilities, since any spin 1/2 density operator 
can be written as 

P = ^ {l + (o^)av + (3^)3^1 , (10) 
and the relation between P (s) and p is linear. 



2 Spin States in the Probability Representation 

We review in this section the probability representation of spin 1/2 states. This 
representation was introduced in refs. [p^ and treats both, pure states and sta- 
tistical mixtures equally via probability distributions. For a general review, see 
ref.0. 

The central point in such a representation is the property of density matrices 
that, in whatever basis they are written, theirs diagonal elements are nonnega- 
tive probabilities. If a spin state is given by a density matrix p'^/^' (written with 
respect to the eigenvectors), one can describe the same state by a rotated 
reference frame density matrix p'^/^^ (u), where u is the set of Euler angles (f>, 
6, V', which defines the rotation (as usual < 4>,ip < 27r, < 9 < tt). 

Let us describe in some detail how such a rotation is done. For the sake 
of simplicity, let us work first with a pure state, described by a spinor ^' = 
X+ 



again with respect to the eigenvectors. In a rotated reference 

X- 

frame, the same state is described by the rotated spinor given by 

(m) = D (m) (11) 

where 



D(u) 



"e i(0-v,)/2 f,™ L-*(<^+'/')/2 



is an irreducible spin 1/2 matrix representation of the rotation group ||T^. In 
components eq. (|l^) reads 

Xs (u)=^i?L/'^ iu)Xm, 
ni 

and it is worthy to note that the matrix elements Z^im^'' (u) are also named 
Wigner functions (for rotations) . 
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It is now straightforward to pass to density matrix. The same state is de- 
scribed by 

p(l/2) ^ ^^t^ 

and rotation gives 

= (u) (* (u))^ = D (m) ^' (D (m) ^')^ (12) 
= D (u) f^-t (D = D (u) (D (u))^ . 

As eq. ( p^ ) is hncar in p^^/^^ , it also applies for mixed states. 

The nonnegative matrix elements PlJ"^^ (u) have the meaning of probability 
distribution of obtaining i when the rotated z-axis spin component is measured 
in this state. By the geometry of Euler angles it is clear that (u) do not 

depend on t/i, so w ~ {9, (f) can be viewed as points in a sphere (equivalently, 
PlJ"^^ (u) does not depend on the rotation itself, but only on the new oriented 
quantization axis, defined by (0, 0)). The central object in this approach is then 
defined 

w(^,u)^p%l^\u), (13) 
which is a positive distribution obeying 

Y,w{%,u)^\, (14) 

i 

for all u (i.e.: marginal probabilities property). It should be viewed as a prob- 
ability distribution for spin 1/2 states. 

2.1 Density Matrix from w{i,u) 

In refs. |p^ w{i,u) probability distributions are worked out for arbitrary spin 
states. With aid of some general properties of SU (2) group, inverse integral 
transformations are obtained, and density matrix is given in terms of w^iju). 
As in the present work we are only interested in spin 1/2 states, we will use a 
more direct tomographic approach. 

As spin Hilbert Spaces are finite dimensional, the complete description of 
physical state can be achieved by using a finite set of real numbers. In partic- 
ular, for two level systems, three real numbers are necessary and sufficient to 
completely describe an arbitrary state. Tomographically, it means that three 
well chosen axis can give the complete knowledge of the state of the system. A 
manifestation of these facts, with a canonical choice of axis, is given by eq. (p^. 



The mean values involved in eq. (10) can be directly obtained by the defini- 
tion of w{i,u) distribution, since 

-u- ( , (15) 



2 / V 2 , 

where Ou = ^ ■ u, and marginal probability property (eq. (p^)). One obtains 

{^)^2w(+^,u)-l. (16) 
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Substitution of eqs. ( p^ in eq. (|l^) then gives 



1 

- I cr„ 



-la. 



(17) 



where we made use of the short notation = w (±i,u(j)), j = x,y,z and 
u{j) the respective </>) pair. An analogous expression holds with wj . In 
matrix form, eq. ( |l7| ) can be written (with aid of eq. ([T^)): 



(1/2) _ 



(18) 



This shows that also w{i,u) contains redundant information, which are 
translated as constraints that a function must obey to describe a quantum 
mechanical state. In this case, the knowledge of three points determine the 
behaviour of the whole function. 

A similar treatment, but which makes use of a nonredundant quorum is made 
by Weigertjl^. In this treatment, the complete state of a spin j is reconstructed 
by measuring 4j (j + 1) (tomographically) independent spin projection mean 
values. 

For completeness, we reproduce here the general formula which gives den- 
sity operator from w {i,u) for any spin value. This formula is deduced in refs. 
from general properties of Wigner's 3j symbols (equivalently, Clebsch-Gordan 
coefficients). It reads 



h=0ma = -j3 



X 



J 

mi 



J ]3 

-mi 



J 



] J3 
-m'2 TO3 



where 



Jl J2 J3 

mi m2 TO3 



dn 



are Wigner's 3j symbols, and 

1 " 



87r2 



d<j) 



smf 



di;. 



2.2 Examples 

In this subsection we compute w {i,u) for the same four examples given in P (s) 
representation. 

i) Spin up in z-axis: In this case, the rotated density matrix reads 



sin f cos f e-»'^ 



■sinf cosfe*'^ 
sin^ f 
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and one gets 



w ( +^,-"1 cos^ ^ = ^ (1 + cos( 



ii) Spin "up" in x-axis: In this case one gets 

w ^+^' ^ ^ ^"'^ ~'~ ^"^^'^°^'^) ' 
w ^"^'""^ ~\^^^ sin^cosi/)) . 

iii) Spin "up" in y-axis: In this case one gets 

=-(1 + sin sin </)) , 

w (^^^'"^ ^ ^ '■^ ^ sin 61 sin (/)) . 

iv) Unpolarized spin: Finally, in this case one gets 

1 \ 1 f I 



— ,u\ — — — w , 
2/2 V 2 



(19) 



This last example simply exhibits the isotropy of the state. The other three 
examples are best understood if one remembers polar spherical coordinates ex- 
pressions. 

In fact, the w distribution for spin 1/2 can be visualized with aid of the 
Bloch ball (i.e. Bloch sphere for pure states and its interior for mixed ones). 
The connection between Bloch vector b and density matrix is made by (see eq. 
®) 

p(^/^) = ll + b-§, (20) 

where b = ^ (^^"^ Euclidean norm lower than or equal to ^. From the 
definition of w one obtains 

«^(±^,") =^±^-^?, (21) 
where u is the unitary vector defined by the polar angles {6, </>). 



3 QUASIPROBABILITY P FROM PROBABIL- 
ITY w 

In this section we explicit formulas which give the quasiprobability distribution 
P (s) (in fact the p{c,b,a) of cq. (|^)) from the knowledge of the tomographic 
probability distribution w {i,u) 
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The composition of eqs. (||) and ( p^ ) directly give 










- 






1 

4' 
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+ wt] 


i 

~ 4' 






^) 


[w^ - 


iWy + 


wt] + 


4' 










+ iWy 


+ wt] 


1 

+ 4' 






^) 


K + 


iWy + 


wj] - 


1 

4' 










-iw+ 


+ w^] 


+ 4' 




J'- 






iui^ + 




4' 




I- 






-iw+ 


+ w~] 


1 

+ 4- 



Those equations are analogous to Radon transform, where Wigner quasiproba- 
bility distribution is written down from tomographic probabilities. 

In this way, the "hidden" quasiprobability P (s) can be obtained directly 
from the knowledge of tomographic probabilities . 

In particular, the four examples worked out previously can be checked, and 
the knowledge of the three Cartesian axis tomographic probabilities allows one 
to completely determine the "joint quasiprobability" P. 

Concluding Remarks 

We have developed a Radon like transformation, in which P quasiprobability 
distribution for spin 1/2 states is written in terms of (some special points of) 
the tomographic probability distribution w. 

The theme of (quasi)probability distributions for spin states is much more 
rich, and besides these two treated here, there are some others "Wigner" distri- 
butions proposed by diverse authors, with different motivations. As an example, 
in ref.|jl^, Chumakov, Klimov and Wolf show that a Wigner quasiprobability 
distribution proposed by AgarwalpOf is a restriction of a Wigner function later 
proposed by WolfpH for generic Lie groups. 

Some peculiar properties of P quasiprobability are evidenced. They are 
related to the operator ordering question, also addressed in ref . pl|] . 

The connection between Bloch vector and w tomographic distribution is also 
made (eq. (l2ll)). It shows a direct geometric interpretation of w in the case of 
two level systems. 
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Appendix: Calculational tools 

In this work we made use of Pauli spin matrices in the form 



1 

1 



-i 

1 



1 
-1 



Conventionally, 



\±x) = 

\±y) = 



1 
1 

71 



{\z)±\-z)}, 
{\z)±z\-z)}, 



and follows 



i+x 1 


+y) {+y 1 


+z) {+z 1 


+x) 


{-x\ 


-y) {-y 1 


+z) i+z 1 


-x) = 


{-x\ 


+y) i+y 1 


-z) {-z 1 


-x) = 




-y) {-y 1 


-z) {-z 1 


+x) = 


{-x\ 


-y) {-y 1 


-z) {-z 1 


-x) = 




+y) i+y 1 


-z) {-z 1 


+x) = 


{+x 1 


-y) {-y 1 


+z) i+z 1 


+x) 


{-x\ 


+y) i+y 1 


+z) i+z 1 


-x) = 


i+x 


+y) i+y 1 


+z) {-z 1 


+x) 


-{-X 


1 -y) {-y 


1 +z) {-z 


\-x) = 


-{-X 


1 +y) i+y 


-z) (+z 


\-x) = 




-y) i-y 1 


-z) {+z 1 


+x) = 




-y) {-y 1 


+z) {-z 1 


+x) 


-{-X 


1 +y) i+y 


1 +z) {-z 


l-x) = 


-{-X 


1 -y) {-y 


-z){+z 


\-x) = 




+y) i+y 1 


-z) {+z\ 


+a;) = 



i(i-i), 
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